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I. INTRODUCTION 


To determine validated sets of elementary chemical reactions for 
use in predictive combustion models, we have chosen to simulate the one 
dimensional, premixed, laminar flame. This approach has the advantage 
that the predicted temperature and species profiles can in principal be 
compared with suitable burner experiments of the same flame. The corre- 
sponding equations are derived by Bird, Stewart and Lightfoot, and by 
Williams,¢ and are given in Appendix A. 


Such a simulation requires as input not only the kinetics informa- 
tion (of our immediate interest), but also thermodynamic and transport 
data. Fortunately the thermodynamics input is by and large well de- 
fined.3,4 In addition, while some transport coefficients are only well 
defined through low temperature (< 1000 K) measurements,° the theory is 
sufficiently developed to allow reasonable estimates to be made at higher 
temperatures.° A theory has been developed for multicomponent 


AH B. Bird, W. S. Stewart and E. N. Lightfoot, Transport Phenomena, 


John Wiley and Sons, NY, (1960). 


ae A. Williams, Combustion Theory, Addison-Wesley, Reading, MA, (1965). 


°p. R. Stull. and H. Prophet, JANAF Thermochemical Tables, 2nd Edvtton, 
NSRDS-NBS-37, June 1971]. 

a Gordon and B. J. McBride, "Computer Program for Calculation of 
Complex Chemical Equiltbrtum Compositions, Rocket Performance, Inet~ 
dent and Reflected Shocks and Chapman-Jouguet Detonattons", NASA-SP- 
273, (1971), (1976 program verston). 


Oy, S. Toulouktan, P. E. Liley, and S. C. Saxena, Thermophystcal 
Properties of Matter, Vol. 3, Thermal Conductivity (Nonmetallte Liquids 
and Gases), IFI/Plenum, NY-Washington, (1970). 


Si 0. Hirshfelder, C. F. Curtiss and R. B. Bird, Molecular Theory of 
Gases ‘and Liquids, 2nd printing, corrected, with notes, Jonn Wiley and 
Sons, NY, (1964). 


mixtures ,6-19 but it is computationally cumbersome. To circumvent this, 
previous workers have generally employed some level of simplification, 11-21 


"CO. 8. Wang Chang, G. E. Uhlenbeck and J. deBoer, Studies in Statistical 
Mechantes, Vol. 2, John Wiley and Sone, WY, (1964). 


The Monchtck, K. S. Yun and E. A. Mason, "Formal Kinetic Theory of 
Transport Phenomena in Polyatomte Gas Mixtures", J. Chem. Phys., 39, 
654-699 (1963). 


L. Monehick, A. N. G. Peretra and BE. A. Mason, "Heat Conductivity in 
Polyatomte and Polar Gases and Gas Mixtures", J. Chem. Phys., 42, 3241- 
5256 (1965), 


L. Monehtck, R. J. Munn and BE. A. Mason, "Thermal Diffusion in Poly- 
atomte Gases: A Generalized Stefan-Maxwell Dtffuston Equation", J. 
Chem. Phys., 45, 3051-3058 (1966). 


G. Dtxon-Lewis, et al., "Flame Structure and Flame Reaction Kinetics", 
Proc. R. Soec., London A 817, 255-263 (1970); A 330,. 199-218 (1972); A 
i 071-584 (1973); and A 346, 261-278 (1975). 
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G. Drxon-Lewts, "Kinette Mechanism, Structure and Properties of Pre- 
mixed Flames tn Hydrogen-Oxygen-Nitrogen Mixtures", Phil. Trans. R. 
Soe., London, 292, 45-99 (1979). 


G. DPsatsaronts, "Prediction of Propagating Laminar Flames tn Methane, 
Oxygen, Nitrogen Mixtures", Combust. and Flame, 88, 217-239 (1978). 


18 


wae Warnatz, "Caleulatton of the Structure of Laminar Flat Flames 1; 


Flame Veloetty of Freely Propagating Ozone Decomposition Flames", Ber. 
Bunsenges Phys. Chem., 82, 193-200 (1978). 


D. B. Spalding and P. L. Stephenson, "Lamtnar Flame Propagation in 
Hydrogen + Bromine Mixtures", Proc. R. Soc. London A, 324, 315-337 
(1971). 


P. L. Stephenson and R. G. Taylor, "Laminar Flame Propagation tn Hydro- 
gen, Oxygen, Nitrogen Mixtures", Combust. and Flame 20, 231-244 (1973). 
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HOR i, Smoot, W. C. Hecker, and G. A. Williams, "Prediction of Propa- 


gating Methane-Atr Flames", Combust. and Flame 26, 323-342 (1976). 


J. Warnatz, "Calculation of the Structure of Laminar Flat Flames ITI; 
Flame Veloetty and Structure of Freely Propagating Hydrogen-Oxygen and 
Hydrogen-Atr Flames", Ber. Bursenges Phys. Chem., 82, 643-649 (1978). 


D. B. Spalding, P. L. Stephenson, and R. G. Taylor, "A Calculation Pro- 
cedure for the Predtctton of Laminar Flame Speeds", Combust. and Flame 
17, 55-64 (1971). 


L. Bledjzan, "Computation of Time-Dependent Laminar Flame Structure", 
Combust. and Flame, 20, 5-17 (1973). 


E. Cramarossa and G. Dixon-Lewis, "Ozone Decomposition in Relation to 
the Problem of the Extstance of Steady-State Flames", Combust. and 
Flame 16, 243-251 (1971). 
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This paper addresses the question: which of the mathematical 
approximations to the multicomponent transport properties provides a 
desirable trade-off between precision and computational effort. Another 
way of phrasing this question is to ask what loss in precision of pre- 
dicted flame speeds and profiles occurs as the mathematical approxima- 


tions to the multicomponent, polyatomic transport expressions are made 
eruder, 


We approach this problem by actually computing the flame properties 
for a given flame with a fixed set of input conditions. The transport 
algorithms are varied and the resulting flame speeds and profiles com- 
pared. The numerical method is discussed in references 22 and 23. We 
have selected the H2-02-Nz system as the candidate flame because it has 
a well characterized set of input parameters. (See Appendix B). 


This paper is divided into the following sections. Section II out- 
lines a sequence of five methods for approximating the transport 
coefficients. Section III shows the numerical results of these approxi- 
mation methods relative to the most accurate method, and Section IV dis- 
cusses their relative precision for a set of five flames. Based on these 
results, we outline in Section V an algorithm which provides a reasonable 
balance between computational effort and theoretical rigor and includes 
a crude approximation for the effects of thermal diffusion. There are 
several appendices that contain sufficient detail to reproduce our com- 
putations. 


II. APPROXIMATION METHODS 


In this section we outline five approximations to the multicom- 
ponent, polyatomic formalism, based on the theory of Wang Chang and 
Uhlenbeck .©7 We start with the most accurate and progressively consider 
cruder approximations. Method I is the three term Sonine approximation 
to the formalism. Dixon-Lewis (DL)24 gives a discussion of this approxi- 
mation, and we will generally follow his approach. Method II involves a 
level of approximation using smaller matrices. Methods III and IV use 


i) Bp) Coffee and J. M. Hetmerl, "A Method for Computing the Flame 


Speed for a Laminar, Premixed, One Dimenstonal Flame", BRL Technical 
Report, ARBRL-~TR-02212, Jan 80. 
ao P. Coffee, "A Computer Code for the Solutton of the Equations 
Governing a Laminar, Premixed, One Dimensional Flame", BRL Memorandum 
Report, to be published. 
aac Dixon-Lewis, "Flame Structure and Flame Reaction Kinetics IT. 
Transport Phenomena in Multteomponent Systems", Proc. Roy. Soc. A 807, 
177-185 (1968). 
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non-matrix techniques. Method V discusses the case of constant trans- 
port parameters. In each case we need to determine the diffusion velo- 
cities Vi, i= 1, 2, ...N, and the heat flux Ges 


A. Method I 


Following (DL), we can write 


N N T 
oT RTD; 0X4 
= _ a ci ot Ee 
Ve 2 py Yi hy ro ax ~ .°| MAY dX (1) 
1=1 1=l1 il 
and 
1. yj Xj Di! 3 
ey, os wg pee Pes 
“i X, © os, 23 oe “Gye ae SP De (2a) 
1 j=l j 1 
or 
V. = V. + W., | (2b) 
iL: 1 Ag 


a notation we will later find convenient. The diffusion coefficients 
can be written in terms of the Sonine expansion coefficients, 


N 
= MX 

k*k o.. 

a 16T = k=l ijk 
Pay “a Gap ° C300 (3) 

j 
and 
T 8MiXq 1 

ms R *i00 ° (4) 


N . a » 
_, WoT 1ji oXj, 8T 1 3 
a 25p ,*, (C300 3x - Bp i00 3x nT) - (9) 
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Similarly, we can write 


=- 4% X. a. - 4 Z X.a ; (6) 


The diffusion velocities and the heat flux are now defined in terms 


1 iba : 
of the 4590? 410° 01? and C599" If we use three terms in the poly- 


nominal expansion, these are defined implicitly by the equations 


(L) a! 1 Pl it ai 1 7 
100 °°? "Noo? “i129 “ee AnTde F161 = “Nol 
4 0.x XX i (7) 
are: Oniky gee: MosMar oara thy 
and 
(1) (chk Ink hk lhk  1hk _lhky T 
igo “2 “noe? “lao * “mio® “lol ~~ “Nol 
= (81, - 8), «++ Sy - ype 0 ++. 0, 0... 0) Ask = 1... N, (8) 


The elements of the 3N by 3N matrix (L) are given in Appendix C. The 
procedures we have used for computing these elements are discussed in 
Appendix D. 


From equation (8), it appears that we must solve N? sets of 
equations to obtain the ee In fact, it has been sewonmentiea oe” 
that only the a terms be computed using this level of approximation, and 
that the c terms be computed using a one term expansion (see Method ie 


However, this is not necessary if we observe that we are only 
interested in the Vj's, and not directly in the c terms. We can convert 
the c-vector in equation (8) to a vector whose elements are in fact the 

16T dX, 


V. by multiplying each equation of (8) by 2Ep 7x «and then summing over 
h. After some algebra we find 

T _ ,16T 3X] 16T oXy T 

(L) (V, aus Vis err eKs ) = (55 3x “71 BSp 8x” O =. Oy 3. ax. 8) 


LS 


Equation (9) shows that we can solve directly for the required V; and 
that the right hand side is not a function of either h or k. 


In solving for the transport in our code, we decompose (L) as the 
product of an upper matrix multiplying a lower matrix. Equations (7) 
and (9) can then be solved efficiently by back substitution. It is com- 
putationally efficient to do this rather than generate another, smaller 
matrix for the Vj using a one term approximation and solving this new 
System. To solve these systems of algebraic equations we have used an 
assembly language algorithm2> which has been found to execute on the BRL 
CDC 7600 about five times faster than the equivalent FORTRAN coding. 


B. Method II 


The above formalism is quite complicated to work with, and so 
further simplifications are almost invariably made. 


For diffusion, we can simplify by taking only one term in the 
Sonine polynomial expansions. This results in 


00,00, . lhk Lhk,T _ 


T 7 
(L ) (Clog bbe CNoOO? = (S471 Ke nee San Onk? h,k = 1... N, (10) 


where the elements of eee are defined in Appendix C. This can be 
rearranged® to give the Stefan-Maxwell equations, namely, 


i ON xayy 
= : ., 0; - Ve), PT ode Ns (il) 
This set of equations is not independent, and the constraint 


N 
: Y.V, = 0 (12) 
eal 


must be used in place of one of the equations in (11). Then the dif- 
fusion velocities can be found by solving a set of N equations in N 


unknowns. This approach gives zero for the thermal diffusion VeLOCLty, 
Wi. 


oan C. Hindmarsh, L. J. Sloan, K. W. Fong, and G. H. Rodrique, "DEC/ 


SOL: Solutton of Dense Systems of Linear Algebrate Equattons", 
Lawrence Livermore Laboratory, UCID-80137, 1976. 
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Thermal conductivity cannot be simplified this easily. However, 
ignoring internal energy, we can obtain the matrix form 


10,10 T 


il tat 
L (A140 a a7 9) = (X,; eer Xi) ; (15) 


: 1 26 


where }. = A =a 2 3 Karn. This form will be valid for a 
fe) : 1 110 


7, EF 
mixture of monatomic gases. 


To define the heat conductivity for a mixture of polyatomic gases, 
we adopt Hirshfelder's* Eucken-type relation,27,5 


poly _ mon, 5 Ai : (14) 
m1xX MIX sip 14+ 2 ii XJ 

am Doe 

jf#i 13 i 


The quantity do in method I is not exactly identical to the usual 
thermal conductivity ee but can be interpreted as the thermal con- 


ductivity of a mixture in which the diffusion forces vanish. 


If the thermal diffusion coefficients are zero, as they are at this 


level of approximation, then A_ = oy. 
O mix 


C. Method IIil 


By making additional assumptions, the Stefan-Maxwell equations (11) 
can be further simplified. A common assumption is that all but the jth 
species move with the same velocity V. Then we find that 


x 

The detatled analysis of Monehtck, et Mt.5? shows that the Di7 should 

be Di,int,j. From a practical point of view we have taken Vz int, j 

= tj 

“Og Muckenfuss and C. F. Curttss, "Thermal Conductivity of Multvtcom- 
ponent Gas Mixtures", J. Chem. Phys. 29, 1273-1277 (1958). 

a 0. Hirshfelder, "Heat Conductivity in Polyatomic, Electronically 
Excited, or Chemical Reacting Mixtures. III", Sixth _Internattonal 
Combustion Symposiun, Reinhold Publishing Corporation, NY, 351-366 (1907). 


TS 


wt ; (15) 


V= = a (16) 


which when substituted into (15) yields the formula recommended by 
Hirshfelder and Curtiss¢8 


non de aD OXG 
“4 = y = XJ Ox’ al 
1 j#i i 


Unfortunately, the expression in (17) does not in general satisfy 
equation (12). One technique to satisfy this constraint is due to Boris 
and Oran.29 They note that if a set of diffusion velocities Vj satisfy 
the Stefan-Maxwell equations (11), then so does the set (Vj + Vc), where 
Ve 1s some constant. The value of Ve is chosen such that the constraint 
(12) is satisfied.* 


The heat conductivity formula employed at this level of approximation 
is taken from Mason and Saxena's30,5 simplification of (14), specifi- 
cally 


al 0. Hirshelder and C. FP. Curtiss, "Theory of Propagatton of Flames 


Part I: General Equattons", Third International Combustion Sympostum, 
Willtams and Wilkins Co., Baltimore, 121-127 (1949). 


oom S. Oran and J. P. Borts, "Detatled Modeling of Combustion Systems", 


to appear tn Progress tn Energy and Combustton Setence. 


A 
Borts and Oran also discuss htgher order approxtmattons to the Stefan- 
Maxwell equations that will not be discussed here. 


°0m A. Mason and S. C. Saxena, "Approximate Formula for the Thermal 


Conductivity of Gas Mixtures", Phys. Flutds, 1, 361-869 (1958). 


Gs 


N 1a 
O i=] l + , : d. LJ oe X;) 


2 i 
j#i * 
where 
—— whiten te oe ike 
bo, = PO + Mayme py RIMIy’2 May’? (19) 
Lj ge = ae Mt 


D. Method IV 

Equation (17) or some analogous form has often been used to compute 
diffusion. However, the usual procedure has been to use (17) only to 
compute V], ..., VN-1. Then Vy is computed from (12). This is less 
accurate than the Boris and Oran procedure, expecially for Vue 


Also, an empirical formula for the thermal conductivity, 


v 


(20) 
1s often used.°!?” 


Method IV is comprised of these common approximations. 


E. Method V 


In the case of a binary mixture the Stefan-Maxwell equations (11) 
reduce to Fick's law. Specifically we have, 


y.v. = -p.. 2% 


a 12 ox (21) 


op, H. Burgoyne and F. Wetnber, "A Method of Analysts of a Plane 


Combustton Wave", Fourth Symposium on Combustion, Williams and Wilkins 
Co., Baltimore, 294-302 (1983). 
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A generalization of (21) can be aaa, and ylelds 


OY; 
a 
a Dam ox? ey 
where 
D. = alae i : (23) 
im » 
jfi 0; 


In the Lennard-Jones formalism Djj «= pie dy_gQ Gl where gf1,1)" 
1s approximately proportional to q-0.17, 32 Since p « T-! it is not un- 
reasonable to assume that p¢ Dij is approximately independent of temper- 
ature. Generalizing, it is often assumed that (p“Dj,,) 1S constant. | 
Likewise for a monatomic gas i « TO.S/Q2,2)*, where Q(452)* « T-0.16 94 
and it is not unreasonable to assume that (pA) is approximately indepen- 
dent of temperature. We now outline a procedure that permits an 
a prtort selection of these quantities. 


For a given flame we know Ty and Yjy, the temperature and the mass 
fractions of the unburned mixture. We also have a chemical kinetics 
scheme and a method of computing the specific heats cpj and specific 
enthalpies hj. Since enthalpy is conserved, i.e., 


, (24) 


we can compute the adiabatic flame temperature Tg. The numerical pro- 
cedure is to guess a trial burned temperature Tt and, using an ODE pack- 
age, find the corresponding equilibrium mass fractions Yjt. We then 
compute the enthalpy of the burned mixture and compare with the unburned 
mixture enthalpy. Tt is iteratively adjusted until the burned and un- 
burned mixture enthalpies agree to within a pre-determined error toler- 
ance. We then accept these values of Tt and Yjt as Tp and Yip. 


Ms al Westenberg, "Present Status of Informatton on Transport Pro- 


perties Applicable to Combustion Research", Combustion and Flame, 1, 
646-358 (1957). 
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At this level of approximation we would like to use an appropriate 


constant (or "global"') Cys. So we will assume that Chi = Cy = constant, 
i = 1... N. Then the relation h, = ht + 4 c . dT becomes h. = h: + 
a 1 ue pl 1 1 
es Cs - By and the mixture enthalpy is given by 
N N 
> Vike = 2 Gh PG. = FD (25) 
By | Sle - : p O 
i=1 1=1 
Substituting (25) into (24) we find 
N 
© 
a NG ts Yip? 
= ‘ (26) 
tT, - ae 


As a heuristic rule, we select T = 0.5 (Tg + Ty) and Yq = 0.5 
(Yyp + Yyy) and then evaluate Dim and pA by using equations (23) and 
(20), respectively. Then p“Dim, 1 = 1,2,..-N-1, and pA are evaluated. 
The diffusion velocity Vy is found from equation (12). 


Others!9-4! have employed the concept of constant a Der, oA and Ch» 
but the method of evaluation varies from author to author. 


III. NUMERICAL RESULTS 


Table 1 shows a summary of the five methods of computing the trans- 
port properties used in this paper. Five H7-07-N2 flames were selected 
and their initial conditions listed in Table 2. The total pressure is 
fixed at one atmosphere for all flames. The computed flame speed for 
each flame as a function of transport method is tabulated in Table 3. 
(The flame speeds for flame A are not corrected to 291 K, as has been 
done.!! If this were done, the value AI for example would be 12.2 cm/s 
instead of 14.1 cm/s.) The values of the flame speeds span a large 
range and for a given flame are essentially independent of the transport 
method. The largest difference between Method I, the most complete for- 
mulation of the transport, and any other method is 16%. (Compare Methods 
I and III, flame D). Note that even Method V gives results that are quite 
close to the much more complex Method I. 
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TABLE 1. SUMMARY OF TRANSPORT METHODS 


Method Remarks 
I 5 terms of Sonine polynomial expansion (3.N by 
3 N matrix); only method that has non-zero thermal 
diffusion. 
II For diffusion, 1 term of Sonine polynomial expansion 


(N by N matrix), Eq. (11) and (12). For thermal 


conductivity, Hirshfelder-Eucken method (N by N 
matrix); Ba: (13) and (14). 


III Diffusion velocities computed from simplified 
Stefan-Maxwell relation, Eq. (17). Each 15 adjusted 
by a common factor, V_, so as to satisfy © Y.V. = 0. 


Thermal conductivity from Mason and Saxena, Eq. (18) 
and (19). 


IV Diffusion velocities computed from the simplified 
Stefan-Maxwell relation, Eq. (17), for N-1 species. 
VN 1S computed from Eq. (12). Empirical thermal 
conductivity formula, Eq. (20). 


2 
V 0 Dim = constant; pA = constant; Cp = constant; 
constants are determined a prtorv. 


TABLE 2. INITIAL TEMPERATURE AND MOLE FRACTIONS 
FOR THE FIVE FLAMES STUDIED 


Flame XH2 X02 XN2 Zu 
A « LSS 0460 5 Oe 7 596 
B . 2000 . 1680 £6320 298 
C . 9000 . 1050 «83350 298 
D . 9000 ~ L000 0 298 
E 6000 4000 0 298 
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Flame 


|r 


14. 


98 


Za2 


378 


892 


i 


TABLE 3. 


375, 


922 


EAN 


FLAME SPEEDS CALCULATED USING 
THE FIVE TRANSPORT MODELS 


Iv v 
14.9 1620 
103 96 
308 291 
402 348 
969 847 


Reproduction of flame speeds is a necessary but not sufficient con- 
dition to judge the relative effectiveness of the transport methods. We 
must also examine the species and temperature profiles of these flames. 
As examples we consider two sets of profiles that exhibit differences 
among the five methods that are as large as any observed. Figures 1 and 
2 show the OH profiles for flame D. Figures 3 and 4 show the H> profiles 
for flame C. As can be seen these profiles are very similar. Almost all 
the other species profiles and all temperature profiles show a greater 
degree of similarity among the five profiles than the examples given. 


In general, our results agree with those reported by Warnatz. ?° The 
exception is flame A, for which Warnatz reported a flame speed of 10.4 
cm/sec, referenced to 291K. Our computed flame speed is about twenty 
percent higher. Based on our results, we infer that Warnatz's use of a 
transport method slightly different from any of those reported here should 
not make that large a difference. We do not know the cause of the 
discrepancy. 


As a check of our code, we computed the flame speed for flame A, 
using the kinetics and transport parameters of Dixon-Lewisl? and Method I. 
The resulting flame speed was 9.2 cm/sec, referenced to 291K. This agrees 
with the computed result reported by Dixon-Lewis and we conclude that our 
code is functioning properly. 


We are comparing the model results produced as we vary the transport 
algorithm and not comparing the model result with experiment. Thus, the 
fact that the use of Warnatz' parameters produces a slightly inaccurate 
flame speed for flame A is of no concern. 


IV. DISCUSSION 


The numerical results shown in Table 3 demonstrate that reliable 
results can be obtained for the H7-07-N2 system even for the case of our 
a prtort determined constant transport method. Note that we cannot infer 
that transport is unimportant! The computed profiles and flame speeds 
can be sensitive to the choice of transport parameters selected. For 
the relative tests of the transport methods here, we have employed the 
same set of species viscosities, thermal conductivities and binary dif- 
fusion coefficients in all cases. We have demonstrated that the method 
used to generate the multicomponent, polyatomic transport coefficients 
is not critical for the H2-07-N> flame. And since this flame is reason- 
ably complex we infer that this result has a high probability of being 
valid for other flames. 


Indeed, we conclude that gross errors detected in comparing the 
results of different models are more likely to be traceable to differences 
in input data rather than to the method of approximating multicomponent 
polyatomic transport properties. 


We cannot completely explain the surprising agreement among the 
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various methods. In this section we will show that our results are con- 
sistent with the few previous studies of transport algorithms. We will 
also discuss some additional numerical results that give some indication 
of the processes involved. 


Warnatz 4 essentially compared Methods II and III for an ozone 
flame and found the differences negligible. However, it was not clear 
that his observations could be generalized to more complex flames. 


In another paper’® Warnatz included thermal diffusion in a simula- 
tion of an Hj-07-N2 flame. He used essentially Method IV together with 
a simple binary formula for thermal diffusion valid only for heavy 
isotopes. © He found that the effect of thermal diffusion so defined was 
to lower the flame speed in a rich flame by about 5%. Species and tem- 
perature profiles were only slightly changed. He did not discuss the 
accuracy of this approximation for a mixture. 


As noted earlier, pene considered the effect of thermal 
diffusion. However, he assumed a particular mass flow and simplified 
chemistry. He found that the thermal diffusion flux of H2 was quite 
large and inferred that effects due to thermal diffusion would be impor- 
tant when considering hydrogen atom diffusion in flame systems. This 
inference is contrary to our findings. Our results also show that, for 
H and Hj, the thermal diffusion velocity Wj can be quite large compared 
to the species diffusion velocity Vj. However, we find that these dif- 
ferences cause only small changes in both the flame speeds and the pro- 
files. This finding supports that of Dixon-Lewis!2 who in a recent paper 
reported that the inclusion of thermal diffusion for near stoichiometric 
H,-air flames lowers the burning velocity by 5 to 6%. 

In an attempt to gain additional insight, we have written an auxilary 
code to directly compare the transport properties. The input required 
for this code are the values of the temperature, the species concentra- 
tions, and the temperature and species gradients. Then the thermal 
conductivity do and the diffusion velocities Vj; are computed using the 
subroutines written for Methods I through IV, and comparisons among the 
values for X45 and Vy are made. By doing this we avoid the complicated 
interplay between transport and chemistry that occurs in flames. 


The auxiliary code was executed using several hundred sets of input 
values that were typical of those observed in the flames reported above. 
These input values were divided into four groups for analysis. Group 1 
consisted of small species and temperature gradients, Group 2 of large 
species and temperature gradients, Group 3 of large species gradients 
and a zero temperature gradient, and Group 4 of small species gradients 
and large temperature gradients. By a large gradient we mean a value 
typical of a flame front. 


The values of thermal conductivity computed with the auxiliary code 
showed little variation, regardless of the input group or approximation 


Ze 


method used. Specifically, the differences between Methods I and II, I 
and III, and I and IV averaged 2, 3 and 5%, respectively. In general, 
the thermal conductivities of the individual species are not this well 
known either theoretically or experimentally. So despite the widely 
different levels of complexity, the four methods produced equivalent 
values for Ao: 


The results of the comparisons for diffusion velocities are more 
difficult to summarize. Differences varied from less than 1% up to 
several hundred percent, depending on the particular input values and 
the particular species considered. Differences of 30 to 40% were common 
and in general were larger for the major species than for the minor species. 


We could not establish why there was a high degree of similarity 
in the flame speeds and profiles even though these were large differences 
in the V. in the auxiliary code. We speculate that part of the reason 


for this is due to the constraint Yr Y.V; = 0, which is satisfied by 


all of our methods. Another constraint on the diffusion velocities in 
the flame code is the fixed chemistry employed. How the non-linear feed- 
back between the chemistry and the transport affects the Similarity in 
the flame speeds and profiles is unknown. 


Nevertheless, some interesting conclusions can be. drawn from these 
comparisons. For instance, it is often assumed that thermal diffusion 
is negligible. This assumption corresponds to that made in Group 3. 
(If the temperature gradient is zero, the thermal diffusion velocities 
W will also be Zero). 


For this group, we find that the values of the diffusion velocities 
computed by Methods I and II agree to about 1%. In contrast, the values 
of the diffusion velocities computed by Methods III and IV differ from 
those computed by Method I by much more (typical differences are around 
40%). This comparison implies that Method II (i.e., the Stefan-Maxwell 
equations) is much more accurate in determining the values for the dif- 
fusion velocities than the non-matrix methods. Note that this conclu- 
Sion is valid only for the case of zero temperature gradient. 


For the other three groups of input to our auxiliary code, we find 
that the diffusion velocities computed by Method III are usually in 
better agreement with those computed by Method I than those computed by 
Method II! This situation can occur because in Method III the errors in 
ignoring thermal diffusion and in approximating the molecular diffusion 
velocities are generally of opposite sign. We conclude that it is a 
waste of effort to use a matrix method (such as the Stefan-Maxwell 
equations) to compute diffusion velocities unless thermal diffusion is 
also considered. 


Method III is generally more accurate than Method IV. Method IV 
can be extremely inaccurate in computing the diffusion velocity for the 
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Nth species. In the cases considered we have found that the values 
computed by Method IV for No differ, on the average, more than 100% from 
the values found by Method I. This is due to the procedure used to 
obtain the diffusion velocity for the Nth species. This procedure tends 
to associate cumulative errors caused by neglect of the constraint, Eq. 
(12), with the Nth species, here No. We have found that Method III (the 
Boris and Oran procedure) is at least as accurate and in some cases 
substantially more accurate than the traditional Method IV for computing 
the V;. Furthermore the Boris and Oran procedure involves very little 
additional computational effort. 


V. A NEW TRANSPORT METHOD 


For flames more complicated than H7-02-Nz we require a transport 
algorithm that is computationally efficient and relatively precise. The 
most exact procedure considered, Method I, can become prohibitively 
expensive for a large number of species. Method V is computationally 
efficient, but we feel caution is necessary in using such a simplified 
model. There does not seem to be a great deal to choose from among 
Methods II, III, and IV and so we have assembled a new method, Method VI. 


As with Methods II-V we require expressions for A and for the 
Vi(= Vy + Wy). For the thermal conductivity we use the simplest formula 
(20), since the exact choice does not appear to be important. For the 
molecular diffusion velocities Vj we use the expression in equation (17). 


For problems with light species and steep temperature gradients, 
the neglect of thermal diffusion is often as important as the differences 
between the computational methods. So below we generate a technique that 
approximates the thermal diffusion contribution to the diffusion velocity, 
W;. We shall first derive an expression for W; of a binary mixture and 
then generalize the results. 


For a binary mixture equation (2a) can be written 


eh 
i XX, 1.2 0x py) dX 


from which we identify 


“= * Dy dent . (28) 
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Since,» from Eqs (LZ). 


. ~ Ul-Y1) OX] | 
Vi. XX Di dx (29) 


— X1X9_ Di : (50) 


we can write, after some algebra, 


gAL 


W = Ri Vy (anT) / (31) 


] 


Theoretical expressions for k)7 have been derived and even in the 
first approximation the expressions are quite complicated.© For the 
special case of heavy isotopes, however, the first approximation simpli- 
fies too 


— IS (ZA* + 5) (Ge* = 5) (My =e 
12 2A*(16A* - 12B* + 55) (M, + M. 


) 
y~ XX, (32) 


Fortunately Eq. (32) reproduces H5-N2 thermal diffusion ratios to within 
30%, and so we use it as a simple but reasonable approximation for Ry>. 
Equation (31) can now be evaluated (in the code). 


Generalizing Eq. (31) we define 


ake 
W. 5 Re VW. = (ent) / sob j (33) 


If an expression for Rjm can be obtained, then Wj of Eq. (33) can be 
evaluated. Chapman and Cowling’’ have derived an approximation for R. > 


a. Chapman and T.G. Cowling, The Mathematical Theory of Non-Untform 
Gases, Third Editton, Cambridge University Press (1970). 
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specifically 


kh _f hk, (34) 


where Rs is given by an expression analogous to Eq. (32). 


From Eq. (32) we can see that the influence of Wj will be the 
greater the larger the mass differences. Normally the thermal diffusion 
ratio does not exceed 0.1 and for the H2-05-Nz system, we have computed 
W; only for the species H and Hj. The neonT tang diffusion vel@eities 


v. = V. + W. do not satisfy the constraint 2 Y.V. = 0. We use the 
1 1 1 oie 


Boris and Oran procedure discussed in Section II to obtain this condition. 


Computationally, the quantity 15(2A* + 5)/(2A*(16A* - 12B* + 55)) 
can be fitted very accurately by a function of the form: a exp (c/T*). 
The quantity (6C* - 5) varies over the range zero to one and could not 
be fitted accurately by a simple expression. The fit for C* given in 
Appendix D is used instead. 


This new method for approximating multispecies transport has been 
exercised in our test flames. Table 4 compares the flame speeds for 
Methods I and VI. The results are so close that we infer that some of 
the errors made in the approximations involved in Method VI are cancelling. 


Figure 5 compares the OH profiles for flame D and figure 6 compares 
the H, profiles for flame C. In general, for the minor species, the 
accuracy of Method VI is comparable to Methods II or III. However, in 
almost every profile involving a major species, the accuracy was improved. 


VI. ADDENDUM 


A few additional comparisons have recently been made. In the spirit 
of approximation exemplified in the discussion of Method V, the further 
assumption is sometimes made that all the Lewis numbers (Le; = A/ cp oDim) 
are equal to unity. Given values for pi and Cp, this assumption defines 
the diffusion coefficients and removes the requirement of supplying 
independent information for them. We tested how well this additional level 
of simplification fared by applying it in the code using flame C. The 
computed flame speed was about 40% lower than for Method I. While the 
major species profiles were still reasonably close to those of Method I, 
they were noticeably less accurate than Method V. The minor species 
profiles showed large differences: the peak concentrations of 0, HO02, and 
H707 were lower than those of Method I by about a factor of two. 


bul 


In an effort to substantiate our earlier inference concerning the 
transfer of these results to other (hydrocarbon) flames we performed some 
checks on the stoichiometric methane-air flame. For this flame we used 
the input parameters listed by Tsatsaronis,*” the kinetics scheme was 
used as listed, and we made no attempt to critically evaluate it. The 
flame speeds computed for Methods I, V, and VI are in the ratio 1.00 to 
0.88 to 0.94, respectively. The profiles for Method VI are more accurate 
than those produced by Method V. These results follow the same trend as 
in the H., -0.-N. flames reported above. 
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TABLE 4. FLAME SPEEDS CALCULATED USING TWO TRANSPORT MODELS 


Flame L VE 
A 14,1 14.0 
B 98 97 
C 292 292 
D 378 402 
B 892 894 


So 
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APPENDIX A. THE FLAME EQUATIONS 


We are interested in the equations that describe a one-dimensional, 
laminar, premixed flame that propagates in an unbounded ideal gas. The 
effects of radiation, viscosity, and body forces are ignored. Since the 
burning velocity is small compared with the local speed of sound, the 
pressure is taken to be constant. The resulting equations (see for 
example. references 2 of text) age: 


Overall Continity: 


ot xR et) 

Continuity of Species: 
p His pu Si = - 2 (py,V,) + RM, i= 1,...N . (A-2) 
and Conservation of Energy: 
aT aT. .B&, 9 (oYiV4) 

oe, aa © puc., = a ax ~ = h, piety nh | : (A-3) 
The boundary conditions are the following (t > 0). For x = -~; 

Rel @my,=V,, Glatell) ¢ (A-4) 
and for x = © ; 

ioe & svi - 0 (i=1,...N) . (A-5) 


The transport properties that we are interested in are the diffusion 
velocities Ve and the heat flux q. 
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APPENDIX B. H,-0,-N, INPUT PARAMETERS 
The required input data for the H2-02-N2 flames are for the most 
part taken from Warnatz (see reference 18 of text). 


The kinetic scheme consists of thirty reactions and has been taken 
in toto from Warnatz' Table 3, which for the most part was taken from 
the Leeds review.! Our interest here is a study of the effects of using 
different transport algorithms, and not a critical evaluation of the 
kinetics. So we have uncritically used Warnatz' mechanism. This 
mechanism is in reasonable agreement with experimental data. 


The enthalpies and heat capacities are evaluated using the sixth 
order polynomial fits of Gordon and McBride (see reference 4 oF text). 
These reproduce the JANNAF (see reference 3 of text) values to within a 
few parts per thousand over the temperature range 300K-3000K. 


The transport parameters are given in Table B-2. These are from 
Warnatz except for §;; for OH, which is from Dixon-Lewis (see reference 
12 of text). The values for o and e/k are obtained from viscosity 
measurements. 


The manipulation of these parameters is discussed in Appendix D. 


Bi L. Baulch, D. D. Drysdale, D. G. Horne, and A. C. Lloyd, Evaluated 


Kinette Data for High Temperature Reactions, Vol. 1, Butterworths, 
London 1972. 
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TABLE B-1. REACTIONS IN THE H,-O,-N. SYSTEM 


REACTION 


OH + H, + H,O + H 


2 


H + HO + OH + H, 


H+ 0, > OH + 0 


O + QH>H + 0, 


O + H., > OH + H 


H + QOH>+ QO + HL, 


OH #4 OH = HO + O 


O + HO > OH + OH 


H +H + M? > H, + M' 


H + OH + M! + HO se 


H+ 0, + M' + HO, + M' 


H + HO. > OH + OH 


H + HO, + H, + 0, 


H + HO, + H,0 + 0 


OH + HO., > HO + 0, 


O + HO., > OH + 0, 


OH + OH +M+ HO, + M 


H + HO. > HO + OH 


HO. “r HO, > HO. + 0, 


H, + MW > hae H+ Me 


O+ 04 HO, + u 


HO, + M" > H +0, + M! 


OH + OH + HO., + H 
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Z 
A* 


Z 


Z 


B 


0 


0 


NO REACT ION 

24 HO, 4: H, > H,0, + H 

25 H + H,0, > H, - HO, 

26 H,0, + M~> OH + OH + M 
27 O + H,0, > OH + HO, 

28 O + HO, > H,0 + 0, 

29 OH + H,0, > HO + HO, 

30 HO, + H,0 > H,0. + OH 


[M] = total concentration 


[M'] = 6.0 [H,O] + 1.0 [H,] + 0.4 [0,] + 0.4 IN 


TABLE B-1. 


REACTIONS IN THE H,-0 


2 
Ax 

7 eBid 

L,~ZE1Z 

i... 2617 

2:8E13 

2.8E13 

1.0E13 


2.8E13 


B 


0 


-N, SYSTEM (Cont'd) 


C 


- 9400 


a 
A ts tn untts of om’ /mole-sec or em® /mole“-see, A= AT” etp (C/T). 


aA 
2.2E13 = 


2.26 & 10°, 
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TABLE B-2, 


Species ofA] 
0 Aud d 
0. 3.46 
H 2405 
H, 2882 
OH ZnlS 
H,0 ZOO 
N. Sez 

HO 


92 10, Ltke O6 


MOLECULAR PARAMETERS USED FOR 
THE DETERMINATION OF TRANSPORT PROPERTIES 


EK] uy [Debye] 
80 O 
107 0 
145 0 
58 0 
80 O 
372 1.844 
97.5 O 
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a [A>] 


APPENDIX C. DEFINITION OF THE MATRIX L 


The matrix L can be partitioned into nine N by N matrices 


00290 20,10 00,01 
fs pi »90 119,10 pi291 (C-1) 
101,00 01,10 pot,01 
where 
pis0d =. 6. - (C-2) 


00,00 _ 16 T [X4x; Mj XpXj| ; 
“47.0 25 P EF Pre M.D. | : ear al, (S-2) 


* 
.1¢ 8 BO2 69 - Ba «) : ; 
Bad ~ SP WM, + M,) oe a (he4) 


13 
peers cn gy QP » (C-5) 
ii So ane 
j#i 
ptOp 08 aq OPatl (C-6) 
ij Hat 
Ot OO o: { POPE = G (a) 
ij ij 
2 
10,10  - 16M; X 10 Gi ot \ . 16T X3Xbh 
ie a” = BBS (; + ar RE) eee | X 
11 ls n 37 R baa 25P bd (M. + Mp)? D-, 
Te, 2. Ol a 2 2% * 5 
= My * Ge Mp wee, ae iM MAS e( * oT 
wha BOR. » che tot 2 | (C-8) 
ik ik ; 
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10,01 
ii 


L, 


10,01 
1j 


L 


Lot 20 
ij 


L 
1) 


and 


L 
aT 


* 
«52 T°M) A1]_ xX Xi. Cl, rot 
2 pT Cs (M, - Mo Ds 6. 
= L 


01,10 _ 


01,01 _ 


_ _16Mz Xi* Cy. rot 


i, int pai My * ME) Vin bap 


J pane ji 


10,01 
ee : 


ART , ex 
PCr aint |e “i, int,k 


| 
12 My A ik XiXk C3. rot 

+2 Gri DL. ¢ 5 ® 

bi h Pir “i, int Sik 
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(C-9) 


(C-10) 


(C-11) 


(C-12) 


(C-13) 


(C-14) 


When solving equations (7) and (6) for the thermal conductivity, 
Ao will turn out to have the units cm2-atm-s~!-K-+. To convert to the 
more usual units of cal-cm-!-s-1!-K-1, we multiply by R/Rg = 1.9872/82.05. 
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APPENDIX D. MANIPULATION OF TRANSPORT DATA 


The quantities required to compute the elements of the matrix L 
used in Method I can be calculated from the transport parameters (Table 
B-2) and the heat capacities (see reference 4 of text). We shall discuss 
the procedures in this appendix. 


The viscosity nj and the binary diffusion coefficients Djj can be 
obtained from Lennard-Jones (non-polar) or Stockmayer (polar) parameters 
ok and ex/k. For polar molecules (here only H20) we also need the dipole 
moment ui. For the binary diffusion coefficient of a polar and a non- 
polar species, we also need the polarizibility a4 of the non-polar species. 
Then nq and Vij are calculated using tabulated values of the collision 
integralsl,2 (see also reference 6 of text). 


* * * 

The quantities A ij» B ij. and C jj do not have exactly their stan- 
dard meanings (see reference 9 of text). However, they can be closely 
approximated by the usual quotients of collision integrals. For 
Lennard-Jones parameters, these are tabulated in reference 6 of text. 
For simplicity, we will assume that at most one of the species is polar 
(in our case H20). Since a polar molecule interacting with a non-polar 
molecule can be represented by a Lennard-Jones potential, all the 
quantities A" ij, B*ij, and C*jj can be cast into the Lennard-Jones for- 
malism. They are functions only of T"4j = T/(eij/R), where ejj is 
obtained from ej and ej by standard combining rules! (see also reference 
6 of text). 


The collision numbers €jj for rotational relaxation of species 1 on 
colliding with species j are difficult to obtain, particularly if i#fj. 
Fortunately, these values are not critical. For simplicity, we will 
assume that €j4 = 47. Parker” has calculated values for £ii for homo- 
nuclear diatomic molecules, as well as the temperature relation 


-l 


co rad la ne 1 
Gaia = Gas 2 | 27, ye + (3 + ‘a | ° (D-1) 


We use his values for 07 and N2. For other types of molecules, we must 
look for experimental values. In this paper, we take these other &jj's 
to be constant. The values for jj are taken from Warnatz (see reference 
14 of text) except the value of £44 for OH (see reference 12 of text). 


Dlp ¢. Retd and 1. K. Sherwood, The Properties of Gases and Liquids, 
2nd Edition, McGraw-Hill, NY, (1966). 


Dé7 Monehick and EB. A. Mason, "Transport Properttes of Polar Gases", 


J. Chem. Phys., 33, 1676-1697 (1961). 
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The quantities Di int, j refer to the diffusion of internal energy, 
which is approximated by Dij. More accurate values are very difficult 
to calculate, and there is no known way of determining them experi- 
mentally (see reference 9 of text). 


To obtain the Ci, rot and Ci,int, we need to properly partition 
Ci,v, the heat capacity at constant volume. Ci,v can be calculated from 
the Gordon and McBride polynominal fits (see reference 4 of cext). Tt 
can be partitioned as 


ie ~ Oi tr * Si eit? : a 


where Ci,tr = 1.5 R. The internal energy can be further partitioned 
into rotational and vibrational energy. Ci,rot = 0 for a monatomic 
molecule, Cj rot = R for a linear molecule, and Ci,rot = 1.5 R for a 
polyatomic molecule. 


For Method I, we do not need to directly compute the separate species 
thermal conductivities 3. These values are required for the other 
methods, so we will discuss them here. 


For monatomic gases, 


(D-3) 


where a hard sphere model has been assumedl (see also reference 6 of 
text). 


For polyatomic gases, there are many procedures. We will follow 
the analysis of Mason and Monchick,1,4 and write 


Hit Cr aeactr ine 7 5 enit 


MUiog c.. +f, C (D-4) 
‘ | 


where 


a G. Parker, "Rotational and Vibrational Relaxation in Dtatomte 


Gases", Phys. Flutds, 2, 449-462 (1959). 
as 2 A. Mason and L. Monehick, "Heat Conductivity of Polyatomte and 
Polar Gases", J. Chem. Phys., 36, 1622-1639 (1962). 


20 


5 10 2 Mi oDiint,i (<3 rot 
= — —— es EE ee ES ee 
Fag ts [ (1 : } wba tt 


37 ns R é,, (rot) 
z (c- i] 
R é... (vib) : (D-5) 
aid. 
fe — cl EPI SInE. 617 & > (1 + “22 —— 
ny NG i > uF ( i ere. Ag CSO) 
OE . 
of 2 ‘s (D-6) 
Oy Sie iw) 


Normally, &4i(vib) >> &4i1(rot), so we will ignore the vibrational 
term. We then refer to &jj(rot) just as —jj. As usual, we approximate 
Di,int,i by the self diffusion coefficient Dii. 


Using either Lennard-Jones or Stockmayer collision integrals, 


Mi oDii 


art 


a ts 
= 1.1999 A yn 1.32. (D-7) 


Equation (D-4) then simplifies to 


AME 1.32c. + 1.77R - .886 Cierot (D-8) 
i ee Si 


As before, the exact value of baa is not critical in the computation. 


For actual use in the code, a fit to the nj, Aj, and Dj; is made of 
the form aTD over an appropriate temperature range, using the logarithmic 
method with proper transformation of weights.» The resulting fits are 
accurate within two percent. 


Don 3. Cvetanovic and D. L. Singleton, Internat. J. Chem. Kinet., 9, 


481-488, 1977; R. J. Cvetanovte and D. L. Stngleton, Internat. J. 
Chem. Kinet., 1007-1009, 1977. 


od 


The quantities A*, B*, and C* are fit over a range 1.0 < T* < 400. 
This covers the range likely to occur in a flame. If the fit is extended 
to values of T* < 1, the accuracy suffers. A* can be fit very accurately 
by a function of the form: a(T*)b. To fit B* and C* accurately, we 
need to use the functional form: a exp (c/T*) (see Table D-1). 


Given the above information, we have developed software to write 


the actual transport subroutine used in the code. 


* * * 
TABLE D-1. FITS TO THE FUNCTIONS A go git 
(a (T*)b or a exp (c/T*)) 


Function a b € 

A. 1.0885 9 .0557E-3 0 

B 1.0786 0 8 .5001E-2 
Cc 95426 0 ~ 14357 
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APPENDIX E. PROGRAM LISTINGS 


We have included here the FORTRAN programs for writing the transport 
subroutines, as well as examples of the subroutines themselves. Some of 
these routines are very long for the H2-02-N2 system (over 3500 lines for 
Method I). So as an example we have considered a binary mixture of H?2 
and N2. This shows the main features of the programs while conserving 
space. 


A copy of the actual flame code used is given in reference 23 of the 
text. Method V (constant transport) is included in this code as the 
default. The appropriate constant transport values for a given flame 
must be read in, but the subroutine itself is not rewritten. 


For our other transport methods, the subroutine is written specifi- 
cally for a given set of chemical species. The procedure is straight- 
forward. We attach to TAPE15 the coefficients for computing the thermo- 
dynamic quantities. We attach to TAPE9 the least squares fits for the 
viscositites, thermal conductivities, and binary diffusion coefficients. 
There are computed in another program using the methods discussed in 
this article. The program also requires the appropriate Lennard-Jones 
(Stockmayer) parameters, the pressure, and for some methods the least 
squares fits for A*, B*, and C*. The flame code actually works with non- 
dimensional forms of the flame equations, so appropriate normalization 
terms for the space coordinate, time coordinate, and temperature must 
also be supplied. 


The transport statements are written onto TAPE3. This tape is read 
and attached to a standard set of initial statements. The corresponding 
FORTRAN subroutine F can then be attached to the flame code. 


Given the mass fractions, the non-dimensional temperature, and 
their first and second space derivatives, F computes the corresponding 
time derivatives. The needed chemistry terms are found by calling the 
Subroutine RT. This subroutine is also written by a FORTRAN program 
(see reference 22 of text). 


For details on the actual form of the equations used and the con- 
nections between F and the rest of the flame code, see reference 23 of 
the text. 


The listing that follows gives the FORTRAN program (MAIN) and the 
corresponding transport subroutines for a binary mixture of H?2 and N2 


(F) for Methods I, II, III, IV and VI. 
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LIST OF SYMBOLS 


‘ (2,2)* ,(1,1)* : Reese 
A ij aT ee » quotient of collision integrals. 

* : ; * * * 
B 13 = - ~ aay?) yy? » quotient of collision integrals. 
‘a OP ee Pee ; , 
Cpe B Rex ja. , quotient of collison integrals. 

iJ 1) L J 

—— = internal heat capacity, eaienoie="_¢,, 

Laie 

a rotational heat capacity, cal-mole >-K~?. 
C. 4. = translational heat capacity, ca'lenolle SK te 
C. a heat capacity at constant volume, sclenole =k. 
ae specific heat, cal-gm t-K7?, 

N 
a; = - rk specific heat of the mixture, cal-gm *-K7?. 
i=l 

= binary diffusion coefficient, ene" 

i = diffusion coefficient of species i into a mixture, en’ -2"* 
ve = multicomponent diffusion coefficient, ened 
D.! = multicomponent thermal diffusion coefficient, en-s"t-cem™>. 
U. news = binary diffusion coefficient for internal energy, fins 
he = specific enthalpy, cal-gm"? 

-16 -] -] 
k = Boltzmann constant = 1.38054X10 erg-K ~-molecule 
Rs = thermal diffusion ratio. 
R. = thermal diffusion ratio between species i and the rest of the mixture. 
M. = molecular weight, em-mole ~~ 
p = pressure, atm. 
q = total energy flux relative to the fluid velocity, wilech “SS 
Rs = rate of production of ith species by chemical reaction, mole-cm™-s +, 
R = gas constant = 1.9872 seine 4K 
_ _ 3 -l ,-1 

R. = gas constant = 82.05 cm -atm-mole “-K -. 


1 


LIST OF SYMBOLS (Cont'd) 


-l 


t = time, s. 
T = temperature, K. 
i. = temperature of unburned mixture, K. 
tT, = temperature of burned mixture, K. 
T, = T/(e,/R). 
* 
Ts = T/Ce, ,/k). 
u = fluid velocity, cm-s7 +. 
Ve = diffusion velocity, cm-s”* 
V. = diffusion velocity due to species gradients, cm-s". 
W = AaPEuStiol velocity due to the temperature gradient, cm-s 
X = spatial coordinate, cm. 
X. = mole fraction, 
Y. = mass fraction, 
ae = mass fraction of species i in the unburned mixture. 
YB = mass fraction of species i in the burned mixture. 
: a 0 
oe i polarizibility, A. 


e./k = Lennard-Jones or Stockmayer parameter, K. 


n, = viscosity, poise = gm-cm +-s~ 

u,; = dipole moment, debye = io" gn gre a" 

Oo, = Lennard-Jones or Stockmayer parameter, he 

ro = thermal conductivity of the mixtures, neglecting diffusion effects, 
cal-cm-l-s-1-xK-1, 

ao = thermal conductivity due to translational energy, eal~on *~s+-xK 

o_o = thermal conductivity due to internal energy, éalecnnas a. 


LIST OF SYMBOLS (Cont'd) 


fe) aT -l 
i = thermal conductivity of a mixture of monatomic gases, cal-cm - 
s-1-K-1, 
poly — ; -l - -1 
Daa = thermal conductivity of a mixture, caul-cm -s -K -. 
o = fluid density, gm-em~. 


bij = rotational collision numbers. 
* * V\* * 
g (sl) : q (2.4) : a heed ; (1,3) = collision integral. 
1) 1J Lj 13 
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